Inducing energy gaps in graphene monolayer and bilayer 
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In this paper we study the formation of energy gaps in the spectrum of graphene and its bilayer 
when both tliese materials are covered with water and ammonia molecules. The energy gaps obtained 
are within the range 20-30 meV, values compatible to those found in experimental studies of graphene 
bilayer. We further show that the binding energies are large enough for the adsorption of the 
molecules to be maintained even at room temperature. 

PACS numbers: 71.15.-m,71.15.Mb,71.20.Tx,73.22.-f, 
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I. INTRODUCTION 

Graphene is a two dimensional system made solely of 
Carbon atoms arranged in a honeycomb lattice. Since its 
isolationii^ in the late 2004, it has attracted a great deal 
of attention due to its unconventional physical properties. 
Together with this wonderful system it was also possible 
to produce its layered relatives, like the graphene bilayer, 
trilayer, etc. 

The physical properties of single layer graphene are 
governed by the elementary excitations around the Fermi 
surface, which in graphene is made of only six points in 
the Brillouin zone, located at momenta values denomi- 
nated K and K' points^. These elementary excitation 
are massless Dirac particles, described by the 2-1-1 Dirac 
equation. The unusual properties of graphene are there- 
fore a consequence of the fact that the emergent physical 
theory of this system is not the usual Schrodinger equa- 
tion but instead a relativistic problem. For the layered 
systems the situation is also equally interesting. For ex- 
ample for the bilayer case, although the low energy spec- 
trum is parabolic, the effective theory is that of massive 
chiral particles^i^ and again not the usual Schrodinger 
equation. 

Both graphene and its bilayer are expected to integrate 
the next generation of nanoelectronic devices^. However, 
both graphene and its bilayer do not present an energy 
gap in the spectrum, preventing large on-off signal ratios 
in possible electronic devices. In the case of graphene, 
there is also a finite conductivity valueii^i^ at the neu- 
trality point, which prevents the pinch-off of the field 
effect transistor. Clearly, what is necessary for having 
full working devices made of graphene is an energy gap 
in the spectrum, as in the usual semiconductors devices. 

There are at the moment two different proposals for 
producing energy gaps in the spectrum of graphene. 
The first is by producing graphene nanoribbons^ i i ^ ^ . 
The theoretical studies^i^^ show that depending on the 



type of terminations a graphene ribbon may have, the 
mechanisms for the opening of a gap will be differ- 
ent. For zig-zag nanoribbons, the existence of edge 
statesi^^ii^ would prevent the existence of an energy 
gap. However, the system finds a way out by inducing 
magnetizatioi*^ at the edges and opening a gap in the 
spectrum. The same happens in a graphene bilayer with 
zig-zag edgesM^. For nanoribbons with armchair edges, 
the gaps are formed due to quantum confinement, and 
the value of the gap depends on the width of the ribbon. 
For ribbons smaller than 10 nm on/off ratios as large as 
10^ have been demonstratcdli. 

Another method of producing a gap in the spectrum 
of graphene is by depositing graphene on top a boron 
nitride (BN)i^. This material is a band gap insulator 
with a boron to nitrogen distance of the order of 1.45 A3, 
(in graphene the carbon-carbon distance is 1.42 A) and 
a gap of the order of 4 eV. It was shown that in the most 
stable configuration, where a carbon is on top of a boron 
and the other carbon in the unit cell is centered above 
a BN ring, the value of the induced gap is of the order 
of 53 meV. Depositing graphene on a metal surface with 
a BN buffer layer leads to n— doped graphene with an 
energy gap of 0.5 eV-^. 

The existence of energy gaps in the spectrum of 
graphene also prevents the occurrence of the Klein 
parado x^^i^" , which would be another limiting source to 
the effective pinch off of the field effect transistor. 

The situation in bilayer graphene is also being unveiled 
experimentally. It has already been demonstrated that 
it is possible to open a gap in the spectrum of the bi- 
layer when this material is deposited on top of silicon 
carbide (SiC) and its exposed face is covered with Potas- 
sium atoms^^. The value of the gap is, in this case, con- 
nected to the amount of doping the potassium atoms 
induce. Using a similar method, it was experimentally 
demonstrated the opening of a gap in a bilayer deposited 
on top of silicon oxide, by covering its exposed face with 
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ammonia molecules^^. In both these methods, the value 
of the gap and the amount of doping are interconnected. 
Recently2^, a device with both a top- and a back- gate se- 
tups became available, which allows to control the value 
of the gap in the bilayer spectrum and the amount of 
charge carrier independently. 

/^From the above experimental studies it is clear that 
covering the bilayer system with some atomic species will 
lead to a gap in the spectrum. These results motivated 
first principle studies^l of gap formation in graphene bi- 
layer leading to calculated gaps in the interval 0.64-3 eV. 

In this paper we study the effect of covering both 
graphene and its bilayer with water and ammonia, show- 
ing that in both cases this covering leads to a gap in the 
spectrum, albeit smaller than those reported in Ref . [13] ■ 
Our results for the induced gaps in the graphene bilayer 
are compatible with those determined for this system by 
measuring the cyclotron masses of the charge carriers^^. 



II. THEORETICAL METHOD 

We are interested in describing the effect that cover- 
ing graphene and its bilayer with water and ammonia 
has in their electronic spectra and in particular in the 
appearance of an energy gap in the spectrum. There is 
strong experimental evidence for the presence of water 
in graphene»2^, and there are working graphene devices 
based on covering one of the surfaces of the graphene 
bilayer with ammonigk^^. Therefore our motivation is to 
understand from first principles how these two types of 
molecules change the electronic spectrum near the Dirac 
point. 

The method we use in our investigation rests upon the 
local density approximation (LDA). In particular, the 
calculations were performed using the density functional 
code AiMPRo2^i22. In what follows we give the relevant 
details of our calculations. 

The Brillouin-zonc (BZ) was sampled for integrations 
according to the scheme proposed by Monkhorst-Pack^^. 
A grid of 8 X 4 X 1 k-points was generated and folded 
according to the symmetry of the BZ. An increase in 
the number of points did not result in a significant total 
energy change. 

We use pseudopotentials to describe the ion cores. 
Lower states (core states) arc accounted for by using the 
dual-space separable pseudopotentials by Hartwigsen, 
Goedecker, and Hutter— . The valence states are ex- 
panded over a set of s-, p-, and d-like Cartesian-Gaussian 
Bloch atom-centered functions, and Fermi- filled using a 
value of fc^T = 0.01 eV. Kohn-Sham states are expressed 
as linear combinations of these basis functions. The ba- 
sis functions were optimized for graphene, water and am- 
monia separately. The calculation of potential terms is 
more efficient in reciprocal space, and therefore the elec- 
tron density is Fourier transformed by using plane waves 
with kinetic energy of up to 300 Ry. 

Graphene was modeled in a slab geometry by includ- 
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FIG. 1: (Color online) Water on the top of graphene: top and 
side views. 

ing a vacuum region in a supercell containing 8 carbon 
atoms. Different sizes were also tested, as well as several 
water or ammonia molecules concentrations per super- 
cell, leading to very different molecular surface densities 
but to only small changes in the gap values. In the normal 
direction (z-direction), the vacuum separating repeating 
slabs has more than 30 A. The size of the supercell in 
the z-direction was optimized to make sure there was no 
interaction between repeating slabs. 

In the calculations, all the atoms were fully relaxed to 
their equilibrium positions. That is to say that the final 
structure of the system is that obtained by energy min- 
imization and is not imposed externally. Consequently 
all the distances found below are a consequence of en- 
ergy minimization of the full structure. 

The electronic band structure for the path K — F — M 
was calculated for all the structures. Since the supercell 
is much larger than the primitive cell, the bands appear 
folded in the graphs. The electronic density of states 
(DoS) was calculated by sampling 12800 points in the 
BZ, and broadened with a 0.1 eV Gaussian width. 

The binding energies were calculated by subtract- 
ing the total energies of the graphene sample and the 
molecule calculated separately from the total energy of 
the graphene plus molecule. 



III. RESULTS 
A. Structure 

Figure [1] shows the relaxed water molecule on the top 
of the graphene. In this calculation, water concentration 
is very low (one water molecule per 8 carbon atoms, or 
a density of about 6 nni^^) and this results in low inter- 
action between the water molecules, meaning that each 
water molecule acts like an independent scattering cen- 
ter. One of the hydrogen atoms binds to a carbon atom, 
and the carbon atoms move from their initial positions, 
specially in the z-direction (as much as 0.05 A). The pres- 
ence of water therefore breaks the A-B symmetry of the 
graphene lattice structure. This, as we shall see, will 
affect the electronic states (see Section Fill Bp . specially 
close to the Dirac point. The same happens when we add 
ammonia. 
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FIG. 2: (Color online) Ammonia on the top of graphene: top 
and side views. 



The ammonia in the top of the graphene (figure [2|) 
has a higher binding energy (see table since it binds 
through two hydrogen atoms, and not through one, as 
the water molecule. Table |T] shows the distances be- 
tween the heaviest atom of the water and the am- 
monia molecules and the closest carbon atom. We 
choose the position of the heavier atoms of the molecules 
as a reference because their position can usually be 
measured using techniques such as low-energy electron 
diffraction (LEED), X-ray photoemission spectroscopy 
(XPS), ultraviolet photoemission spectroscopy (UPS), 
and electron-energy-loss spectroscopy (EELS)^. Our re- 
sults show a much higher binding energy than the ones 
obtained by other authors^. We believe this is because 
those authors impose restrictions to the rotation of the 
molecules on the top of graphene; we observed that our 
molecules never relaxed to those positions, but rotate in- 
stead to more favorable orientations. 

When a second layer of graphene is added, a bilayer is 
obtained. The usual stacking of the graphene bilayer is 
the Bernal staking, although there are situations where 
rotational defects may also be present^. In the Bernal 
stacking two of the carbon atoms, each belonging to dif- 
ferent layers, are positioned exactly on top of each other. 
For these atoms we have found that they move away from 
the plane as much as 0.004 A. The distance between lay- 
ers (measured as the distance between two carbon atoms 
that are superimposed in the Bernal stacking) is 2.6 A. 
This is lower than the interplain distance on graphite. 
We believe this is due to the absence of layers on both 
sides of the bilayer. 

The bilayer has a binding energy per carbon atom of 
0.21 eV, and both the water and the ammonia molecules 
are less bound to the bilayer than to graphene. This is be- 
cause in the bilayer the carbon atoms near the molecules 
also have to share electrons with the underlying carbon 
layer. In fact, the hydrogen of the molecule that is closer 
to the surface binds to a carbon atom that has no other 
carbon atom below. In every case the binding energy is 
high and so desorption by thermal evaporation at room 
temperature is not expected. This is consistent with the 
experimental finding that graphene usually has adsorbed 
water—. 

We also have done calculations for higher molecular 
concentrations. In these cases, there is strong interaction 




FIG. 3: (Color online) Water on the top of bilayer graphene: 
top and side views. 




FIG. 4: (Color online) Ammonia on the top of bilayer 
graphene: top and side views. 



between the molecules that cause their re-orientation, in 
such a way as to align some of the hydrogen atoms to 
the oxygen/nitrogen in order to make hydrogen bonds. 
For water, two molecules per 8 carbon atoms covers the 
graphene surface; one molecule more and a second water 
layer starts to form. For ammonia, two molecules per 8 
carbon atoms is enough for a second layer to be formed. 



B. Electronic structure 

The electronic band structure for the graphene with 
water and ammonia molecules on the surface is shown 
in figure El Since the wc are using a supercell, the Bril- 
louin zone is different from the standard one, and that 
is reflected on the shape of the bands. But the abso- 
lute location of points K and M are the same. In both 
cases, a localized band appears distinctly above the va- 
lence band of the graphene. It corresponds to the HOMO 
level of the molecules, which is clearly below the Fermi 
level. The LUMO of the molecules is very much above 
the Fermi level. This means that the charge transfer in 
the interaction between the molecules and the surface of 



TABLE I: Distances between oxygen and nitrogen and the 
closest carbon atom of graphene 



System 


Atoms 


Distance (A) 


Graphene 


O— C 


2.79 


Graphene 


N— C 


2.91 


Bilayer 


o— c 


2.83 


Bilayer 


N— C 


2.97 
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TABLE II: Binding energies for the adsorbed molecules of 
ammonia and water (per molecule), and for the bilayer (per 
carbon atom). 



System Binding 


energy (eV) 


Graphene + H2O 


1.94 


Graphene + NH3 


3.58 


Bilayer 


0.21 


Bilayer + H2O 


1.30 


Bilayer + NH3 


2.86 





d) 



FIG. 5: (Color online) Electronic band structure for the sys- 
tem a) graphene plus water, b) detail near the Dirac point, c) 
graphene plus ammonia, d) detail near the Dirac point. The 
clean graphene band structure appears as a shadow in the 
background. 



graphene will be small, and the interaction must be of 
dipole or Van-der-Waals type. 

The bands near the Dirac point can be seen in detail 
in figure [5] and have a similar form for both types of 
molecular coverage of the carbon layers. A gap opens 
and it is shifted to one side of the K point. This happens 
because there is a symmetry lost due to the displacement 
of the carbon atoms interacting with the molecules. 

Table [ni] shows the calculated values for the band gap. 
It is well known that local density functional theory, while 
expected to give a realistic account of the ground state 
structure, it only gives an approximate description of un- 
occupied states in a system. In particular the energy gap 
is underestimated, being typically 60% off the real value. 
This should be taken into account when analyzing the 
bands and the gap values in detail. In principle a gap in 
the graphene spectrum is not expected, except when the 
A-B symmetry of the lattice in broken. As our calcula- 
tions show, there is a vertical distortion of the graphene 
lattice which leads precisely to the breaking of the A-B 
symmetry and therefore to the opening of a gap. This 
breaking of the symmetry lattice translates in a spectrum 
that is hyperbolic near the Dirac point, where the bands 
are still linear for finite energies, except close to the Dirac 
point, where they become parabolic in shape. 

The electronic band structure for the bilayer with wa- 
ter and ammonia molecules on the surface are shown in 
figure [SI Since in our calculations we have the two layers 
much closer than in graphite, which is the distance usu- 



ally considered in other calculations, it is natural that we 
have the second band at higher energies than other au- 
thors, due to an increase of the tight-binding parameters 
in the direction perpendicular to the graphene layers. In 
the case of water on the top of the bilayer there is no flat 
band above the bilayer valence band, while there is one 
in the case of ammonia. 

As in graphene, covering the bilayer with cither water 
or ammonia leads to the opening of gap which is shifted 
to one side of the K point, which is consistent with the 
results obtained by Guinea et alS^. These authors found 
that the gap in bilayer graphene does not reach a mini- 
mum at the K point due to a "Mexican hat" dispersion 
at low energies. Although the resolution of our figures 
is not enough to see "Mexican hat", its presence in the 
band structure was verified by the authors. The gaps 
(see Table HIT)) that the molecules generate on the bilayer 
are much larger than the ones on graphene. 

It is known that the unbiased bilayer has parabolic 
bands at the Dirac point, with the valence and the con- 
ducting band touching each other. The presence of the 
molecules opens up a gap in the spectrum, but while in 
the case of water the gap is very near the Dirac point, 
in the case of ammonia the gap is indirect and far from 
the K point. This last result is a direct consequence of 
the type of molecule used and can not be obtained within 
tight binding calculations. 

It is known22, that a difference in potential along the 
z-direction induces a gap in the electronic band structure 
for the bilayer. The water and the ammonia molecules 
have both a electric dipole which will produce an electric 
field on the bilayer. In this case, the difference in electric 
potential between the layers is 0.198 V for the system 
with ammonia and 0.027 V for the system with water. 

The gap we obtained with ammonia is consistent with 
the gap calculated using tight binding and a selfconsis- 
tcnt determination of the electric potential^^. Yet, our 
calculations lead to a electric potential between graphene 
layers interacting with the ammonia molecule of 0.198 V. 
This value, when included in the cquation^^ 
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gives a much larger gap. The discrepancy may be ex- 
plained by the fact that the expression is not valid for 
indirect gaps, and it does not include the distortion of 
the bilayer lattice. The gap for the bilayer with water 
estimated from the electric potential between the layers 
and using Eq. ^ (based on a tight binding model) gives 
consistent results with those determined from the band 
structure calculation. 

In our case, besides the electric potential induced by 
the molecule, there is also a distortion in lattice both 
in the case of graphene and its bilayer. Our calculations 
show that the distortion of the lattice in the bilayer tends 
to counteract the molecule's dipole, but still leading to a 
band gap. 

The electronic density of states (DoS) of the systems 
considered is shown in figure [T] The molecules create a 
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a) 




FIG. 6: (Color online) Electronic band structure for the sys- 
tem a) bilayer plus water, b) detail near the Dirac point, c) 
bilayer plus ammonia, d) detail near the Dirac point. The 
clean bilayer graphene band structure appears as a shadow in 
the background. 



TABLE III: Calculated band gaps for the systems considered. 
The gap shown for the bilayer+NHg is indirect. 



System Band 


gap (meV) 


Craphene + H2O 


18 


Graphene + NH3 


11 


Bilayer + H2O 


30 


Bilayer + NH3 


42 



lot of structure in the density of states, the most relevant 
are the strong peaks that correspond to very flat bands 
localized around the molecule. Except for the water on 
the bilayer, the molecules generate a very high peak on 
the DoS at higher energies than the first peak of the occu- 
pied states for the clean surfaces. The presence of these 
peaks leads to an increase of the absorption at energies 
of a few eV. 



IV. CONCLUSIONS 



The adsorption of water and ammonia on the top of the 
graphene and graphene bilayer was studied. The position 
and orientation of these molecules relative to the surfaces 
was obtained as well as the binding energies. The high 
values of the binding energies (of the order of the eV) 
indicates that desorption of these molecules cannot be 
obtained at room temperature, and a surface coverage is 
expected if the surfaces are exposed to water or ammonia. 



We have showed that a gap opens up in the band 
structure when molecules adhere to the surface, both on 
graphene and on the bilayer, and its value is of the or- 
der of a few tens of meV. The mechanism for opening an 
energy gap may be different for the graphene single and 
bilayer, since the random distribution of the molecules 
on the single layer may not break the A-B symmetry in 
average while the gap-opening in the bilayer is related to 
the symmetry between top and bottom layers. 



We therefore have showed the existence of a feasible 
method to produce graphenes with gaps in their spec- 
trum, which is not restricted to the manipulation of 
nanoribbons. We found that the gap is in general dis- 
placed away from the Dirac point. 



An increase in the absorption coefficient for energies 
above about 2 eV due to the presence of the molecules is 
predicted from the electronic density of states. 





b) 



d) 



FIG. 7: (Color online) Electronic density of states for the 
system a) graphene plus water, b) graphene plus ammonia 
compared with graphene; c) bilayer plus water, d)bilayer plus 
ammonia compared with bilayer 
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